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Chapter 1 

The dark secrets of gaseous nebulae 



1.1 Emission line nebulae 



The existence and distribution of the chemical elements and their isotopes is a consequence 
of nuclear processes that have taken place in the past in the Big Bang and subsequently in 
stars and in the interstellar medium (ISM) where they are still ongoing (Pagel, 1997). A 
large body of our knowledge of the distribution and production of elements in the universe 
rests on observations and analyses of photoionized gaseous nebulae. Ionized and heated 
by strong ultraviolet (UV) radiation fields, photoionized gaseous nebulae glow by emitting 
strong emission lines (Osterbrock and Ferland, 2005). They are therefore also commonly 
named emission line nebulae. 

Examples of emission line nebulae include H ll regions, planetary nebulae (PNe) and the 
broad and narrow emission line regions found in active galactic nuclei (Fig. H ll regions 



are diffuse nebulae found around newly formed young, massive stars and trace the current 
status of the ISM. Giant extragalactic H ll regions, sign posts of massive star formation 
activities, are amongst the most prominent features seen in a gas-rich, star-forming galaxy. 
In some galaxies, the star forming activities are so intense that the whole galaxy becomes a 
giant H ll region. Such galaxies are called H ll or starburst galaxies and are observable to 
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Figure 1.1: Examples of emission line nebulae. From left to right: 

a) (HST images obtained by B . Balick and collaborators; c.f. 

Ihttp:// www.astro.washington.edu/users/balick/WFPC2/index.html); b) The H ll region 
M4 2 (the Urion INebula; Hubble Heritage image obtained by U. K. O'Dell and S. K. Wong, 
c.f. |http: / /hubblesite^org/newscenter/archiye/releases/lQOS /45 ) ; c) NGC604, a giant extra- 
galactic H II region in the outskirts of the Local Group spiral galaxy M 33 (Hubble Heritage 
image obtained by H. Yang, c.f. http://heritage.stsci.edu/2003/30/supplemental.html; d) 
The starburst galaxy IZwl8 (based on HlS'iy Wl''h'(J2 data obtained by Fl. ykillman). The 
linear sizes of these objects differ vastly, ranging from about a tenth of a parsec (Ipc = 
3.262 lyr = 3.086 x 10^^ m) for PNe and M42 to several hundred parsecs for NGC604 and 
IZwlS. 



large cosmic distances. 

PNe are amongst the most beautiful objects in the sky and, arguably, the queen of the 
night. They were given the name by William Herschel (Herschel, 1785) based on their distinct 
structures and, for some of them, nearly circular and overall uniform appearances resembling 
the greenish disk of a planet. They have however nothing to do with a planet (but see later) 
and are in fact expanding gaseous envelopes expelled by low- and intermediate-mass stars 
in late evolutionary stage after the exhaustion of central nuclear fuel at the end of the 
asymptotic giant branch (AGB) phase. Owing to their relatively simple geometric structure, 
a nearly symmetric shell of gas ionized by a single, centrally located white dwarf, PNe are 
ideal cosmic laboratories to study the atomic and radiative processes governing cosmic low 
density plasmas. PNe have played and continue to play a central role in formulating the 
theory of photoionized gaseous nebulae. Representing a pivotal, albeit transient evolutionary 
phase of low- and intermediate-mass stars (the overwhelming majority in a galaxy), PNe play 
a major role in the galactic ecosystem - in the constant enrichment of metals, formation and 
destruction of molecules and dust grains and in the recycling of gas in the ISM. Today, 
studies of PNe have gone far beyond the objects themselves. PNe are widely used to trace 
the kinematics of the host galaxies and the intracluster stellar populations. They have even 
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been successfully utilized to measure the Hubble constant of the cosmic expansion. 

Recent progress in observational techniques, atomic data and high-performance computa- 
tion have enabled reliable measurements and analyses of lines as faint as one millionth of H/3, 
including weak optical recombination lines (ORLs) from abundant heavy elements (C, N, O, 
Ne and Mg) and coUisionally excited lines (CELs) from rare elements, such as fluorine and 
s- and r-process elements. This allows one to address some of the long-standing problems in 
nebular astrophysics as well as opening up new windows of opportunity. 

In this contribution, I will briefly review the development of the theory of photoionized 
gaseous nebulae, highlighting some of the key events. I will then present some recent devel- 
opments of deep spectroscopy of PNe and H ll regions, concentrating on observations of faint 
heavy element ORLs. I will show that there is strong evidence that nebulae contain another 
previously unknown component of cold (about 1,000 K), high-metallicity plasma, probably 
in the form of H-deficient inclusions embedded in the warm (about 10,000 K) diffuse nebula 
of "normal (i.e. near solar ) composition" . This cold gas emits essentially all the observed 
flux of heavy element ORLs, but is too cool to excite any significant optical or ultraviolet 
CELs and thus invisible via the latter. The existence of H-deficient gas in PNe and probably 
also in H II regions, not predicted by the current stellar evolution theory, provides a natural 
solution to the long-standing dichotomy between nebular plasma diagnostics and abundance 
determinations using CELs on the one hand and ORLs on the other, a discrepancy that is 
ubiquitously observed in Galactic and extragalactic PNe as well as H ll regions. 

1.2 The founding of the theory of photoionized gaseous nebulae 

The apphcations of the principle of "Chemical analysis by observations of spectra" ex- 
pounded by G. Kirchhoff and R. W. Bunsen (1860) by W. Huggins and W. A. Millar to 
the analyses of stellar and nebular spectra in 1860s heralded the rise of astrophysics. In 

an accompanying paper to their monumental work on stellar spectra, Huggins and Millar 
presented their first visual spectroscopic observations of eight PNe, including the Cat's Eye 
Nebula NGC6543 in Draco (Huggins and Millar, 1864). Instead of the dark (absorption) 
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Fraunhofer lines observed in the spectra of the Sun and other stars, they found bright emis- 
sion hnes in the nebular spectra and concluded that the nebulae must consist of " enormous 
masses of luminous gas or vapour" . While the Fraunhofer F line (H/3 at 4861 A) did appear 
in emission, the two bright, nearby lines AA4959, 5007 were not Fraunhofer lines at all. The 
availability of dry photographic plates (light-sensitive silver bromide salts held in a gelatine 
emulsion on the glass), replacing the old wet collodion plates, made it possible to record 
long exposures of nebular spectra and more nebular emission lines were revealed in the blue 
and ultraviolet (UV) wavelength regions. In 1882, Huggins successfully photographed a UV 
spectrum of the nearest bright H II region, the Orion Nebula, and detected a strong emission 
line near 3728 A (Huggins, 1882). 

By the late 1920s, dozens of nebular lines had been detected and their wavelengths accu- 
rately measured (Wright, 1918), yet most of them remained unidentified and were attributed 
to some hypothetical element "nebulium" . The big breakthrough in understanding nebular 
spectra came in 1927 from Ira Bowen who, stimulated by a heuristic speculation by H. N. 
Russell that the nebulium lines must be due to abundant "atoms of known kinds shining 
under unfamiliar conditions" such as in gas of very low density (Russell, Dugan and Stewart, 
1927), identified eight of the strongest nebular lines as being due to the forbidden transitions 
from the low excitation meta-stable states of the ground electron configurations of singly 
ionized oxygen, nitrogen and doubly ionized oxygen (Bowen, 1927a, b, c). For the first time, 
the physical processes in gaseous nebulae could be understood. This important discovery 
paved the way for future studies of nebular structure and chemical composition. 

Great strides forward in nebular spectral observations were made starting in 1910 and 
through the 1930s, powered by a number of technological inventions, such as the Schmidt 
camera (Schmidt, 1932), the high-efficiency blazed diffraction grating (Wood, 1910) and the 
image sheer (Bowen, 1938). Deep exposures revealed faint lines from the refractory elements 
such as potassium, calcium, silicon, magnesium and iron, demonstrating that nebulae are 
qualitatively made of similar material as stars (Bowen and Wyse, 1939). In the following 
three decades, the structures and the underlying physics governing photoionized gaseous 
nebulae had been worked out quantitatively, including photoionization and recombination 
(Zanstra, 1926, 1927; Stromgren, 1939), heating and coohng (Zanstra, 1931; Baker, Menzel 
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and Allcr, 1938; AUer, Baker and Menzel, 1939; Spitzer, 1948), recombination and collisional 
excitation (Baker and Menzel, 1938; Menzel, Aller and Hebb, 1941; Seaton, 1954; Burgess, 
1958; Seaton, 1959; Seaton and Osterbrock, 1957), and elemental abundance determinations 
(Bowen and Wyse, 1939; Menzel and Aller, 1941; Wyse 1942; Aller and Menzel, 1945). 
Origin of PNe, as descendants of red giant stars was also understood (Shklovsky, 1956). 
Photoionization models incorporating all known physics were constructed (Hjellming, 1966; 
Goodson, 1967; Harrington, 1968; Rubin 1968; Flower, 1969) and the models reproduced 
observations well. The theory of photoionized gaseous nebulae as it stood in late 1960s 
was nicely summarized in Astrophysics of Gaseous Nebulae by Donald Edward Osterbrock 
(1974). 



1.3 CELs and ORLs — the dichotomy 



While the theory seemed well estabhshed and solid, there were dark clouds hovering on the 
horizon. One concerned the measurement and interpretation of weak nebular emission lines, 
and the other regarded the possible presence of significant temperature inhomogeneities in 
nebulae and their effects on nebular abundance determinations. 

Except for a few lines excited under specific environments, such as the Bowen fluores- 
cence lines (e.g. O III AA3133, 3341, 3444; Bowen, 1934, 1935), strong lines radiated by 
photoionized gaseous nebulae fall into two categories, recombination lines (RLs) and CELs. 
Hydrogen and helium ions, the main baryonic components of an ionized gaseous nebula, 
capture free electrons and recombine, followed by cascades to the ground state. During this 
process, a series of RLs are radiated (e.g. H a A6563, H /? A4861, He i AA4472, 5876, 6678 and 
He II A4686). The ground electron configuration of multi-electron ions of heavy elements 
yields some low excitation energy levels (within a few eV from the ground state, such as 
0++ 2p^ '^Po.1,2) ^^2, ^Sq), which can be excited by impacts of thermal electrons, which have 
typical energies of ~ 1 eV in photoionized gaseous nebulae of solar composition. Follow-up 
de-excitation by spontaneous emission yields the so-called CELs. Quite often, those lines 
are electron dipole forbidden transitions, so they are commonly called forbidden lines (e.g. 
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Figure 1.2: Top-left: Schematic diagram showing the four lowest spectral terms of O^^ 
formed by the 2p^ ground and 2p^ electron configurations. The levels are labeled with ex- 
citation energy Tgx (in Kelvin) and critical density A^c (in cm~^); Bottom- left: Schematic 
Grotrian diagram of O ll, illustrating the O ll recombination line spectrum produced by 
recombination of O"*"^; Top-right: Photoionization and recombination of hydrogen. Pho- 
toionization heats the gas whereas recombination cools the gas; Bottom-right: Emissivities 
of the collisionally excited [O ill] A5007 optical forbidden line and 52/im far-IR fine-structure 
line, and of the recombination lines H/3 A4861 and O ll A4649, as a function of electron tem- 
perature Te for selected electron densities. The curves are labeled with logarithms of electron 
densities. Emissivities of recombination lines depend only weakly on electron density under 
typical (low-density) nebular conditions. 
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[O II] AA3726,3729, [O lll] AA4959,5007, [N ll] AA6548,6584 and [Ne ill] AA3868,3967). Re- 
combination of free electrons to bound states of hydrogen and helium ions also yields weak 
nebular continuum emission. For example, recombination to the H I n = 2 state yields the 
near UV nebular continuum Balmer discontinuity at A < 3646 A. Recombination of heavy 
element ions, followed by cascading to the ground state, also produces RLs. However, for 
typical cosmic composition, i.e. approximately solar, even the most abundant heavy element 
oxygen has a number density less than one thousandth of hydrogen. Thus heavy element 
RLs are much weaker, at the level of a few thousandth of H/3 or less, and are observable 
generally only in the optical. Those lines are therefore often called ORLs. Sample ORLs 
from abundant heavy element ions that have been well studied and discussed in the current 
contribution include C ll A4267, C ill A4649, N ll A4041, O I A7773, O ll A4089, O ill A3265, 
Ne II A4392 and Mg ll A4481. 

Except for IR fine-structure lines arising from ground spectral terms, emissivities of CELs 
have an exponential dependence on electron temperature, Te, e(X+\ A) oc T~^^'^exp{—Ecx/kTe) 
(the Boltzmann factor), where E^x is excitation energy of the upper level of transition emit- 



ting wavelength A by ion X+' following collisional excitation by electron impacts (Fig. 1.2). 
At low densities, « Nc, i.e. electron density A^^e much lower than the upper level's crit- 
ical density Ac (for any energy level above the ground, a critical density A'c can be defined 
such that above which collisional de-excitation dominates over spontaneous radiative decay 
in de-populating the level; c.f. Osterbrock & Ferland, 2005), we have e(X+', A) oc A^(X''"')Ae, 
where A^(X"''') is number density of ion X+'. At high densities, A"e >> Ac, emission of CELs 
are suppressed by collisional de-excitation and e(X+\A) oc A^(X+^). Unlike CELs, emissiv- 
ities of RLs increase with decreasing T^ by a power law. Under typical nebular conditions 
(A'e « 10*^cm-3), e(X+',A) oc T-"A^(X+'+i)A'e, where a ~ 1 and N(X+'+^) is number 
density of the recombining ion X"'''"*'^. 

In a ground breaking work, Arthur B. Wyse published very deep spectra that he obtained 
with I. S. Bowen using the Lick 36-inch reflector for a number of bright PNe (Wyse, 1942). 
About 270 spectral lines were detected. Many were weak permitted transitions from abun- 
dant C, N and O ions. In the Saturn Nebula NGC 7009 alone, "... more than twenty lines 
or blends of O ll have been observed, free of blending with lines of other elements. There 
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are good grounds for the assumption, previously made, that these hues originate in electron 
captures by O ill ions, just as the Balmer lines originate in electron captures by H ll ions, 
and that therefore the relative intensities of the oxygen and hydrogen lines give a measure of 
the rates of recombination, and therefore of the relative abundance, of the two kinds of ions 
.... The important thing is that for this nebula a sufficient number of oxygen lines has been 
observed to eliminate all doubt of the correctness of their identifications, even though they 
are very faint; and also that they permit the direct estimate of the relative abundance of 
H II and O ill ions, by a method that is relatively independent of assumptions regarding 
electron density and velocity distribution or their variations throughout the nebula." 

Prom the O ll ORL strengths, Wyse deduced 0/H abundances which were much higher 
than the values obtained by Menzel and AUer (1941) from analyses of the [O ill] forbidden 
lines. He concluded "The discrepancy, then, between the relative abundance found by Menzel 
and Aller, on the one hand, and from the recombination spectra, on the other, is of the order 
of 50 for NGC 7027 and of 500 for NGC 7009." Shortly after the pubhcation of this prophetic 
article, Wyse was called to serve in the War and died tragically on duty the night of June 8, 
1942 in a disastrous accident over the Atlantic Ocean, off the New Jersey coast, at the age 
of 17 days shy of 33. 

As stressed by Wyse, ionic abundances deduced from intensities of heavy element ORLs 
relative to H/3, a method based on comparing hues of hke to hke, have the advantage that 
they are almost independent of the nebular thermal and density structures. In contrast, 
ionic abundances deduced from the intensity ratio of the coUisionally excited, much stronger 
[O III] AA4959,5007 forbidden lines relative to H^, have an exponential dependence on the 
adopted nebular electron temperature. Theoretically, ionic abundances deduced from heavy 
element ORLs should thus be more reliable, provided the lines can be measured accurately. 
Unfortunately, latter development showed that accurate flux measurements for faint nebular 
emission lines were NOT possible after all with the technique available then, i.e. spec- 
trophotography, due to the non-linearity of photographic plates. Via detailed comparisons 
between the observed fluxes of H I and He ll RLs and continua and those predicted by 
the recombination theory, it became clear that spectrophotographic observations systemat- 
ically overestimated intensities of faint lines, by as much as over a factor of ten (Seaton, 
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1960; Kaler, 1966; Miller, 1971; Miller and Mathews, 1972). That spectrophotographic mea- 
surements of faint lines cannot be trusted seemed to be further supported by work in the 
1980s that contrasted C++/H+ ionic abundances deduced from the coUisionally excited C ill] 
AA1907,1909 intercombination hues and from the faint C ll A4267 ORL (c.f. Barker, 1991 
and references therein), suggesting that the intensity of the faint C ll A4267 hue had either 
not been interpreted correctly or been grossly overestimated (Rola and Stasihska, 1994). 

In another worrisome development. Manual Peimbert showed that if nebulae are non- 
isothermal and have (localized and random) temperature fluctuations, then deduced from 
the [O III] nebular to auroral line intensity ratio (A4959 -|- A5007)/A4363 will overestimate the 
average emission temperature of the AA4959,5007 nebular lines and of H/3 at 4861 A, leading 
to underestimated ionic abundance ratio deduced from the (A4959 -|- A5007)/H/9 

ratio (Peimbert, 1967). He presented evidence pointing to the presence of large temper- 
ature fluctuations in PNe and H II regions by finding that Te's derived from the Balmer 
jump, Tc(BJ), arc systematically lower than those derived from the [O III] forbidden line 
ratio, re([0 III]) (Peimbert 1967; Peimbert, 1971). Given the weakness of nebular contin- 
uum emission, measuring the Balmer jump accurately was no easy task and his results were 
disputed (Barker, 1978). Theoretically, while some systematic spatial temperature varia- 
tions undoubtedly occur within a nebula, due to changes of ionization structure and cooling 
rates as a function of position induced by varying ionization radiation field and density 
distribution, no mechanisms are known capable of generating large, localized temperature 
fiuctuations, certainly nothing of the magnitude implied by Peimbert's measurements, which 
yield a typical value of 0.055 for the temperature fiuctuation parameter, t"^, or fiuctuations 
of an amplitude of 23 per cent. 

The advent of linear, high quantum efficiency, large dynamic range and large format 
CCDs in the 1980s made it possible for the first time to obtain reliable measurements of 
faint emission lines for bright nebulae. Meanwhile, the completion of the Opacity Project 
(Seaton, 1987; Cunto et al., 1993) has allowed the atomic data necessary to analyze those 
spectral features, specifically their effective recombination coefficients, to be calculated with 
high accuracy. Liu and Danziger (1993) obtained CCD measurements of the Balmer jump 
for a sample of PNe and found that re(BJ) is indeed systematically lower than Te([0 ill]) 
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obtained for the same object. They deduced that on average, = 0.035, smaller than 
that found earlier by Peimbert (1971) but still significant, enough to cause the 0"'"^/H"'" 
abundance ratio derived from the AA4959,5007 forbidden lines to be underestimated by a 
factor of two. 

Liu et al. (1995) presented high quality IPCS and CCD optical spectrophotometry for 
the legendary Saturn Nebula NGC 7009 as well as new effective recombination coefficients 
for O II ORLs. Nearly a hundred O ll ORLs were measured, yielding an 0^^/H+ ionic 
abundance that is consistently higher, by a factor of ~ 4.7, than the value deduced from 
the strong [O III] forbidden hues AA4959,5007. The close agreement of results deduced 
from a large number of O ll ORLs from a variety of multiplets of different multiplicities, 
parities and electron configurations, vindicates the reliability of the recombination theory 
and rules out measurement uncertainties or other effects, such as reddening corrections, 
line blending or contamination of ORLs by other excitation mechanisms (e.g. fluorescence 
or charge transfer reactions) as the cause of the large discrepancy between the ORL and 
CEL abundances. Further analysis of the carbon, nitrogen recombination spectra and of 
the neon recombination spectrum (Luo et al., 2001) show that in NGC 7009, abundances 
of these elements derived from ORLs are all higher than the corresponding CEL values by 
approximately a factor of 5. If one defines an abundance discrepancy factor (adf ) as the ratio 
of ionic abundances X'+/H+ deduced from ORLs and from CELs, then in NGC 7009 adf ~ 5 
for all the four abundant second-row elements of the periodic table: C, N, O and Ne. In 
NGC 7009, the Balmer discontinuity of the hydrogen recombination spectrum yields Te(BJ) 
= 8, 100 K, about 2,000 K lower than forbidden line temperature Te([0 ill]) = 10, 100 K (Liu 
et al., 1995). The difference yields a Peimbert's value of 0.04, or temperature fluctuations 
of an amplitude of 20 per cent. 



^Mathis and Liu (1999) analyzed the observed relative intensities of the [O III] AA4959,5007 and the much fainter A4931 
nebular lines and demonstrated that accurate measurements have been achieved over a dynamic range of 10,000. 
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Figure 1.3: Left: 0/H abundances deduced from ORLs plotted against those derived from 
CELs. The diagonal dotted line denotes x = y. The insert is a histogram of adf(0++/H+). 
In the scenario of single composition nebulae (c.f. next section for an alternative interpre- 
tation), uncertainties in 0/H abundances, caused by observational and interpretive errors 
(e.g. atomic data) are expected to be typically less than 0.05 and 0.1 dex, respectively, for 
CEL and ORL results; Middle: Te([0 ill]) versus Te(BJ). Typical uncertainties of Te([0 ill]) 
and Te(BJ) are 5 and 10 per cent, respectively. The diagonal dotted line denotes x = y. 
With re(BJ) = 900 K and Te([0 ill]) = 8820 K, Hf 2-2 falls off the left boundary of the plot; 
Right: logadf(0"'""'"/H"'") plotted against Te([0 ill]) — Te(BJ). The solid line denotes a linear 
fit obtained by Liu et al. (2001b), prior to the discovery of the very large adf in Hf2-2. 
(Adapted from Liu, 2006a) 



1.4 Deep spectroscopic surveys of ORLs 



Is the dichotomy observed in NGC 7009 between ORLs and CELs for nebular plasma diag- 
nostics and abundance determinations ubiquitous amongst emission line nebulae? What is 
the range and distribution of adf's for individual elements? Are there any correlations be- 
tween the adf and other nebular properties (morphology, density, temperature, abundance, 
and age etc.) or properties of the ionizing source? What are the physical causes of the 
dichotomy? To address those issues, several deep ORL spectroscopic surveys of faint heavy 
element ORLs have been carried out. So far over a hundred Galactic PNe have been sur- 
veyed, plus dozens of Galactic and extragalactic H II regions (e.g. Esteban et al., 2002; 
Tsamis et al, 2003a,b, 2004; Liu et al, 2004a,b; Wesson et al., 2005; Wang and Liu, 2007; 
For a complete list of references, please refer to a recent review by Liu, 2006a). Detailed 



comparisons contrasting ORL and CEL analyses show that (Fig. 1.3; c.f. also Liu, 2006a) 



1) Ionic abundances deduced from ORLs are always higher than CEL values, i.e. adf > 1. 
Adf peaks at 0.35 dex but with a tail extending to much higher values. About 20 and 10 
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Figure 1.4: Left: The optical spectrum of Hf 2-2, which shows a record adf(0"'""''/H+) = 71 
and an extremely low re(BJ) of 900 K. In the lower panel, also shown are two synthetic 
recombination line spectra of O ll and N ll, respectively. Right Monochromatic images of 
Hf2-2 in the light of [O ill] A5007 and Ha A6563, respectively. (Adapted from Liu et al., 
2006) 



per cent of nebulae exhibit adf's higher than 5 and 10, respectively. For example, in the 
bright Galactic disk PN NGC6153, adf = 9.2 (Liu et al., 2000), whereas in the bulge PN 
Ml-42, adf = 22 (Liu et al., 2001b). In the most extreme object discovered so far, IIf2-2, 



adf(0++/II+) reaches a record value of 71 (Fig. 1.4); 



2) While adf varies from object to object, for a given nebula, C, N, O and Ne all exhibit 
comparable adf's (thus both GEL and ORL analyses yield compatible abundance ratios, 
such as G/0, N/0 and Ne/0, provided lines of the same type, ORLs or GELs, are used 
for both elements involved in the ratio). Objects showing large adf's also tend to have high 
helium (ORL) abundances. However, magnesium, the only third-row element which has 
been analyzed using an ORL, shows no enhancement, even in high-adf objects (Barlow et 
al, 2003); 

3) Excluding metal-poor nebulae in the Galactic halo and in the Large and Small Mag- 
ellanic Glouds (LMG and SMG, respectively), oxygen abundances deduced from GELs for 
Galactic H ll regions and PNe fall in a narrow range compatible with the solar value. In 
contrast, ORLs yield much higher abundances, more than ten times solar in some cases. 
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4) Similarly, while the [O III] forbidden line ratio yields values of Te in a narrow range 
around 10,000 K, as one expects for a photoionized gaseous nebula of solar composition, the 
Balmer discontinuity yields some very low temperatures, down to below 1,000 K. In fact, the 
discrepancies in temperature and abundance determinations, using ORLs/continua on the 
one hand and CELs on the other, seem to be correlated - objects showing large adf 's also 
exhibit very low T^s; 

5) Large, old PNe of low surface brightness tend to show higher adf's. In addition, 
spatially resolved analyses of a limited number of bright, extended nebulae of large adf's 
show that ORL abundances increase towards the nebular center, leading to higher adf's near 
the center. 



1.5 Evidence of cold, H-deficient inclusions 

What causes the ubiquitous, often alarmingly large, discrepancies between the ORL and 
GEL plasma diagnostics and abundance determinations? Does the dichotomy imply that 
there are fundamental flaws in our understanding of the nebular thermal structure, or that 
we do not even understand basic processes such as the recombination of hydrogenic ions? 

Can it be temperature fluctuations as originally postulated by Peimbert (1967)? This 
conjecture implicitly assumes that the higher abundances yielded by ORLs represent the true 
nebular composition as they are insensitive to temperature and temperature fluctuations. 
The discovery of nebulae exhibiting extreme values of adf and the abnormally high metal 
abundances imphed by the observed strengths of ORLs casts however serious doubt on this 
paradigm. Given that PNe are descendants of low- and intermediate-mass stars, the very 
high ORL oxygen abundances recorded in objects of extreme adf's, if real and representative 
of the whole nebula, are extremely difficult to understand in the current theory of stellar 
evolution and nucleosynthesis. 

Strong evidence against temperature fluctuations as the cause of the dichotomy between 
ORLs and CELs is provided by Infrared Space Observatory (ISO) measurements of mid- and 
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Figure 1.5: Left The far-IR spectrum of NGC 6153; i^z^'/it: Comparisons of ionic abundances 
of NGC 6153 deduced from ORLs, and from UV, optical and IR CELs. (Adapted from Liu 
et al., 2000) 



far-IR fine-structure lines, such as the [Ne ill] 15.5-/im and [O ill] 52- and 88-/im (Liu et 
al., 2001a). Although collisionally excited, IR fine-structure lines have, unlike their optical 



and UV counterparts, excitation energies less than ~ 1,000K (Fig. 1.2), and are therefore 
insensitive to temperature or temperature fluctuations. If temperature fluctuations are in- 
deed at work, leading to an overestimated Te([0 ill]) and consequently an underestimated 
AA4959,5007 abundance, then one expects a higher abundance from the 52- and 

88-/im fine-structure lines comparable to the value yielded by O ll ORLs. ISO measure- 



ments and subsequent analyses reveal, however, otherwise. Fig. L5 shows that in the case of 
NGC 6153, all CELs - UV, optical and IR, regardless of their excitation energy and critical 
density, yield ionic abundances that are consistently a factor of ~ 10 lower than ORLs. 

To account for the multi-waveband observations of NGC 6153, Liu et al. (2000) pos- 
tulated that the nebula contains a previously unknown component of high-metallicity gas, 
presumably in the form of H-deficient inclusions embedded in the diffuse nebula of "nor- 
mal" (i.e. about solar) composition. Due to the efficient cooling of abundant metals, this 
H-deficient gas has an electron temperature of only ~ 1,000K, too low to excite any opti- 
cal or UV CELs (thus invisible via the latter). Yet, the high metallicity combined with a 
very low electron temperature make those H-deficient inclusions powerful emitters of heavy 
element ORLs. In this picture, ORLs and CELs yield discrepant electron temperatures and 
ionic abundances because they probe two different gas components that co-exist in the same 
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Table 1.1: Comparison of Te's deduced from CELs and from ORLs/continua'^ 



IN tiUUicI 




T (\n tttI'i 




ic[ne I) 
(K^ 




NGC 7009 


4.7 


9980 


7200 


5040 


420 


H 1-41 


5.1 


9800 


4500 


2930 


<288 


NGC 2440 


5.4 


16150 


14000 




<288 


Vy 1-2 


6.2 


10400 


6630 


4430 


3450 


IC 4699 


6.2 


11720 


12000 


2460 


<288 


NGC 6439 


6.2 


10360 


9900 


4900 


851 


M 3-33 


6.6 


10380 


5900 


5020 


1465 


M 2-36 


6.9 


8380 


6000 


2790 


520 


IC 2003 


7.3 


12650 


11000 


5600 


<288 


NGC 6153 


9.2 


9120 


6000 


3350 


350 


DdDm 1 


11.8 


12300 


11400 


3500 




Vy 2-2 


11.8 


13910 


9300 


1890 


1260 


NGC 2022 


16.0 


15000 


13200 


15900 


<288 


NGC 40 


17.8 


10600 


7000 


10240 




M 1-42 


22.0 


9220 


4000 


2260 


<288 


NGC 1501 


31.7 


11100 


9400 






Hf 2-2 


71.2 


8820 


900 


940 


630 



" Adapted from Liu (2006a). Te(He l)'s were derived from the He I A5876/A4472 and 
A6678/A4472 line ratios and were typically accurate to 20 per cent (c.f. Zhang et al., 
2005). Values of Te(0 ll) were deduced from the O ll A4089/A4649 line ratio and had 
typical uncertainties of 30 per cent. 



nebula but have vastly different physical and chemical characteristics. Empirical analysis of 
Liu et al. (2000) as well as follow-up ID photoionization modeling (Pequignot et al., 2002) 
shows that a small amount of H-deficient material, about one Jupiter mass, is sufficient to 
account for the strengths of heavy element ORLs observed in NGC 6153. 

The increasingly lower values of Balmer jump temperature Te(BJ) found for nebulae of 
larger and larger adf's: 6,000 K in NGC 6153 (adf = 9.2), 4,000 K in M 1-42 (adf = 22) 
and 900 K in Hf2-2 (adf = 71), provide the smoking gun evidence that nebulae contain 
two regimes of vastly different physical properties. Further evidence is provided by careful 
analyses of the He I and heavy element recombination line spectra, which show that the 
average emission temperatures of the He I and O ll ORLs are even lower than indicated by 
the H I Balmer discontinuity (Liu 2003). In general, it is found that, Te(0 ll) < Te(He l) < 
Te(BJ) < Te([0 III]) (Table[Ll] c.f. Liu, 2006a and references therein), as one expects in the 
dual-abundance scenario proposed by Liu et al. (2000). 
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Detailed 3D photoionization models of NGC6153 with and without H-deficient inclusions 



have been constructed by Yuan et al. [2010, in preparation; c.f. Figs. 1.6 1) and b)] using 
MOCASSIN, a Monte Carlo photoionization code capable of dealing with nebulae of arbitrary 
geometry and composition (Ercolano et al. 2003a). In their models, the main nebula was 
modeled with a chemically homogeneous ellipsoid of "normal composition" (i.e. about solar 
as yielded by the CEL analysis). To mimic the bipolar shape of the nebula, the density 
in the ellipsoid was allowed to decrease from the equator to the poles. In addition, to 
reproduce the strengths of low-ionization lines, such as [C l] AA9824,9850, [N l] AA5198,5200, 
[N II] AA6548,6584, [O l] AA6300,6363 and [O ll] AA3727,3729, an equatorial ring of the same 
chemical composition but of a higher density was added. The presence of a high-density 
torus is supported by a high resolution spectrum obtained with the Manchester Echelle 
Spectrograph mounted on the Anglo-Australian Telescope. The spectrum, centered on Ha 
and the [N ll] AA6548,6584 lines and obtained with a long-slit oriented in PA = 123 deg and 
through the central star, revealed two high velocity emission spots at the positions of the 
bright shell, one on each side of the central star. The spots are particularly bright in [N ll] 
and have respectively blue- and red-shifted velocities relative to the Ha emission from the 
nearby nebular shell. The model fits the CEL measurements well except for a few lines [c.f. 



the top panel of Fig. 1.6 d)]. The model still underestimates the [C l] A9850 line by a factor of 
three even after the introduction of the high-density torus. The model also underestimates, 
by about a factor of two, the fluxes of the N ill] A1744 and [Ar ill] 8.9- and 21.8-/im lines. The 
latter three lines were however only marginally detected, by the lUE and ISO, respectively, 
and the discrepancies may well be due to measurement errors. The largest discrepancy is 
found for the [Ne ll] 12.8-/im line where the model underestimates the observation by more 
than a factor of five. 

The chemically homogeneous model however fails to reproduce the strengths of all heavy 
element ORLs by a factor of ten. To account for them, metal-rich inclusions were added, in 
the form of clumps of H-deficient material. The clumps were distributed spherically sym- 
metrically with a radial number density profile that reproduced the O ll and C ll ORL 
surface brightness distributions deduced from the long-slit observations (Liu et al., 2000). 
Given that the H-deficient inclusions have not been resolved, even with the HST/STIS spec- 
troscopy with a slit width of 0.2 arcsec (Yuan et al., 2010), they must have dimensions smaller 



1.5. EVIDENCE OF COLD, H-DEFICIENT INCLUSIONS 



17 



HST images 




r=9007K 

2.1 
1 ,9 

1.7 #=0.998 
1.5 M= 0.243 M 

1 .3 

11 r=815K 
u.a - -4000 cm- 
0!5 #=0.002 

0.0031 M 



Figure 1.6: a) From top to bottom: HST monochromatic images of NGC6153 in Ha and 
O III] A5007 (north is up and east to the left); monochromatic images predicted by the 
Dest-fit model incorporating H-deficient inclusions; and the model distributions of electron 
density and temperature. 
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Figure 1.6: - continued, b) Comparison of model line intensities with observations for the 
best-fit chemically homogeneous model (top) and model containing H-deficient inclusions 
(bottom). In the bottom panel, blue and red color bars denote, respectively, contributions 
of line fluxes from the cold, H-deficient inclusions and from the hot, diffuse gas of "normal 
composition" . 



1.5. EVIDENCE OF COLD, H-DEFICIENT INCLUSIONS 19 

than ~ 0.2arcsec, or 300 AU at the assumed distance of 1.5 kpc of NGC6153. The smallest 
size of inclusions that can be modeled with the MOCASSIN was limited by the memory 
available per processor of the computer employed. For NGC6153, a grid of 48^ cells, each of 
1/4 arcsec by 1/3 arcsec, was used to model 1/8 of the nebula using a SGI Altix 330 cluster of 
8 dual-core processors and 8 Gb memory per processor. Models that doubled the resolution 
were also constructed using the NASA/Ames Columbia Supercomputer. Results from one 



of these high resolution models are shown in Fig. 1.6 1) and b). The best-fit dual- abundance 
model incorporating H-deficient inclusions matches all available observations reasonably well 
and successfully reproduces strengths of heavy element ORLs. With a total mass of only 
one Jupiter mass, or just over 1% that of the whole nebula, a helium abundance 5 times 
and CNONe abundances 40 - 100 times higher than the main diffuse nebula of "normal 
composition", the H-deficient inclusions account for approximately, 5, 35 and 100% of the 
observed fluxes of H I, He I and heavy element ORLs, respectively, but produce essentially 
nil optical and UV GEL emission. Gooled to about 800 K by IR fine-structure lines of heavy 
elemental ions, notably the [Ne ll] 12.8-/im, [Ne ill] 15.5-/im and [O ill] 52-/im lines, the 
inclusions are important contributors of those lines, so for the H I recombination continuum 
Balmer discontinuity (about 30%). The inclusions even dominate the [Ne ll] 12.8-/im line 
flux. Note that gas densities in the H-deficient inclusions are not well constrained at the 
moment, partly due to the lack of suitable atomic data and diagnostic tools (see below, §1.8). 
This leads to some uncertainties in the deduced total mass of the H-deficient material. In 
the current treatment, the inclusions have densities and temperatures that are roughly in 
pressure equilibrium with the surrounding diffuse medium of higher temperatures and lower 
densities. Physically, one envisions individual inclusions to be optically thick such that they 
can survive long enough to have observational effects. Proper modeling of optically thick 
knots will however require a separate set of fine grid for each of them to resolve the ionization 
and thermal structures, and thus demand even more computing resources. 

The average elemental abundances for the whole nebula of NGG6153, including the dif- 
fuse nebula as well as H-deficient inclusions embedded in it, are 0.102, 3.9 x 10~^, 4.4 x 10~^, 
7.3 X 10~^ and 2.4 x 10"'* for He, G, N, O and Ne, respectively. Except for helium, elemen- 
tal abundances of the "normal" component are close to those deduced from the empirical 
method based on GELs (Liu et al., 2000), and are about 30 per cent lower than the average 
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abundances for the whole nebula. For helium, the abundance of the "normal" component 
is about 40 per cent lower than derived from the empirical method based on the observed 
strengths of helium ORLs, and actually comes closer to the average helium abundance of the 
whole nebula. That NGC 6153 has a lower helium abundance than implied by the empirical 
analysis is supported by a measurement of the near IR He I A10830 line, for which excitation 
from the 2s meta-stable level by electron impacts dominates the emission (Yuan et al., 
2010). 

Helium lines observable in the optical and UV are all dominated by recombination ex- 
citation, and their strengths can be strongly enhanced by a small amount of H-deficient, 
ultra-cold plasma posited in the nebula. If H-deficient inclusions are indeed fully responsible 
for the ORL versus GEL dichotomy ubiquitously found in PNe and H ll regions, then their 
presence may have a profound consequence on the helium abundances deduced for those 
objects, as well as on the primordial helium abundance determined by extrapolating the 
helium abundances of metal-poor H ll galaxies to zero metallicity - the inferred primordial 
helium abundances from standard analyses would be too high (Zhang et al., 2004). 

Corroborative evidence that ORLs arise from distinct regions of very low temperature 
plasma is provided by analyses of emission line profiles. Spectroscopy at a resolution of 
6kms~^ of NGC 7009 shows that O ll ORLs are significantly narrower than the [O ill] for- 
bidden lines (Liu, 2006b). Observations at an even higher resolution of 2 kms~^ of NGC 7009 
and NGC 6153 yield complex multi-velocity-components and differing line profiles for the O ll 
ORLs and for the [O ill] forbidden lines. The profiles of the [O ill] A5007 nebular and A4363 
auroral lines also differ significantly. These results suggest that ionized regions of vastly 
different temperatures (and radial and thermal velocities) co-exist in those nebulae (Barlow 
et al., 2006). Zhang (2008) shows that the differences in the observed line profiles of the O II 
ORLs and the [O ill] CELs are hard to explain in a chemically homogeneous nebula. 
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Figure 1.7: Left: The optical images of AbellSO, showing the H-deficient knots near the 
nebular center (adapted from Borkowski et al., 1993); Right: The optical spectrum of the 
polar knots J 1 and J 3 showing prominent ORLs of C, N, O and Ne ions (adapted from 
Wesson et al., 2003). Note the faintness of H/3 at 4861 A. 



1.6 Spectroscopy of PNe harboring H-deficient inclusions 



H-deficient clumps were previously known to exist in a rare class of old PNe, including 
Abell30, 58 and 78, IRAS 15154-5258 and 18333-2357 (e.g. Harrington, 1996). They are 
identified as PNe that experience a last helium shell flash that brings them back to the AGB 
to repeat the PN evolution stage, the so-called "born-again" PNe (Iben et al., 1983). In 
Abell 30, HST imaging reveals a host of knots embedded near the center of a round, limb- 



brightened faint nebula of angular diameter 127arcsec (Borkowski et al., 1993; c.f. Fig. 1.7). 
They include two point-symmetric polar knots along PA = 331° at angular distances 6.66 
and 7.44 arcsec from the central star, plus a number of knots loosely delineating an equatorial 
disk or ring. Deep optical spectra of knots J 1 and J 3 were obtained by Wesson et al. (2003). 
The spectra show prominent recombination lines from CNONe ions, remarkably similar to 



those seen in other PNe of large adf's, such as Hf2-2 in Fig. L4, except that in Abell 30 the 
strengths of those CNONe ORLs relative to H/3 are much higher. Detailed ORL analyses 
show that the ORLs are indeed emitted under very low temperatures of only a few hundred 
Kevin. The result is corroborated by detailed 3D photoionization modeling (Ercolano et al., 
2003b) as well as by a similar analysis of the H-deficient knot found at the center of Abell 58 
by Wesson et al. (2008). 
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1.7 Origins of H-deficient inclusions 



The origin of H-deficient material is not well understood and their existence in PNe, and 
possibly also in H II regions (see, for example, Esteban et al. 2002 and Tsamis et al. 2003a), 
is not predicted by the current theory of stellar evolution. For PNe harboring H-deficient 
central stars, such as the Wolf-Rayet and PG 1159 stars, while the scenario of a single post- 
AGB star experiencing a last helium shell flash seems to be able to match the photospheric 
abundances of the star, not all PNe exhibiting large adf 's have an H-deficient central star 
(e.g. NGC7009). The process and mechanism by which H-deficient knots might be ejected 
in such systems are not understood. More importantly, detailed ORL abundance analyses 
by Wesson et al. (2003, 2008) show that the H-deficient knots found in AbeU30 and Abell58 
are oxygen-rich, not carbon-rich as one expects in the scenario of "born-again" PNe. 

The precise coUinearity of the two polar knots in Abell 30 with the central star, to within 
5arcmin, is hard to explain by single star evolution, and suggests instead the action of 
a bipolar jet from, for example, an accretion disk in a binary system (Harrington, 1996; 
De Marco, 2008). In this respect, it is interesting to note that Abell 58 is known to have 
experienced a nova-like outburst in 1917 (Clayton and De Marco, 1997). Hf2-2 has also 
been found to be a close binary system with a period of 0.398 571 d (Lutz et al., 1998). By 
comparing the known properties of Abell 58 with those of Abell 30, Sakurai's Object and 
several novae and nova remnants. Wesson et al. (2008) argue that the elemental abundances 
in the H-deficient knots of Abell 30 and Abell 58 have more in common with neon novae than 
with Sakurai's Object, which is believed to have undergone a final helium fiash. 

An alternative to the scenario of H-deficient inclusions being ejecta of nucleo-processed 
material is that they derive from metal-rich planetary material, such as icy planetesimals, 
that once orbited the PN progenitor star. As the star enters the PN phase by ejecting the 
envelope and evolves to become a luminous, hot white dwarf, the strong stellar winds and 
UV radiation fields begin to photoionize and strip gas off the planetesimals, now embedded 
in the PN (Liu, 2003, 2006a). As analyses show that only a few Jupiter masses metal-rich 
material is required to explain the observed strengths of ORLs, the idea may not be so 
eccentric as it first looks. Evaporating proto-planetary disks around newly formed stars 
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(such as the proplyds found in the Orion Nebula) may hkewise provide a natural solution to 
the problem of ORL versus GEL dichotomy similarly found in H ll regions. 



1.8 The need for new atomic data 



Atomic data relevant for the study of emission line nebulae, including collision strengths 
and effective recombination coefficients, are generally calculated for a temperature range 
5,000 - 20,000 K, typical for photoionized nebulae of solar composition. The finding that 
heavy element ORLs emitted by the H-deficient knots in the "born-again" PNe Abell 30 and 
Abell 58 arise from ultra-cold plasma of only a few hundred Kevin, as well as the compelling 
evidence accumulated so far that points to the presence of similar cold H-deficient inclusions 
in other PNe, call for new atomic data, in particular the effective recombination coefficients 
applicable in such low temperature environments. New plasma diagnostics that use only 
ORLs, not CELs, need to be developed so that physical and chemical properties of the H- 
deficient knots/inclusions, such as temperature, density, elemental abundances and mass, 
can be determined reliably, a prerequisite to unravel their origins. 

The ratios of the intensities of ORLs from states of different orbital angular momenta show 
some temperature dependence and thus can be used to measure the average temperature 
under which the lines are emitted (Liu, 2003). In addition, as pointed out by Liu (2003), 
the dependence on density of the relative populations of the fine-structure levels of the 
ground spectral term of the recombining ions, such as O^"*" 2p^'^Po,i,2, leads to variations of 
the relative intensities of recombination lines as a function of density. In the case of O ll, 
for example, lines from high-/ states, such as the 3p'^Dy2 ^ 3s^P5/2 A4649 transition (c.f. 



Fig. 1.2) weakens as density decreases due to the underpopulation of the O^^ 2p^^P2 level 
under low densities. This effect opens up the possibility of determining electron density 
using ORLs. 



New effective recombination coefficients, calculated down to temperatures of ~ 100 K and 
taking into account the dependence on density of level populations of the ground states of the 
recombining ion, have now been carried out for the O ll (Bastin and Storey, 2006) and N ll 
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Figure 1.8: Lejt: Schematic diagram illustrating the low temperature di-electronic recombi- 
nation of O II via the low-lying fine-structure auto-ionizing states. Note that for ni > 67 and 
77,2 > 37, the doubly-excited O ll 2p^(^Pi)ni/i and 2p^(^P2)n2/2 states fall above the O ill 
ground level 2p^ ^Pq. Similarly, for > 47, O ll 2p^('^P2)n2^2 states have energies above the 
O III first excited fine-structure level 2p^'^Pi. For temperatures of a few hundred Kelvin, 
di-electronic recombination via those low-lying fine-structure auto-ionizing states becomes 
an important process; Right: Loci of the O II recombination line ratios A4076/A4089 and 
A4076/A4069 for different T^s and A'^e's (based on data provided by Dr. P. J. Storey). 



(Fang, Storey and Liu, in preparation) recombination spectra. Fig. 1.8 plots loci of the O ll 
/(A4076)//(A4089) and /(A4076)//(A4070) intensity ratios for different temperatures and 
densities. The diagram also illustrates the importance of di-electronic recombination via the 
low-lying fine-structure auto-ionizing states under very low temperatures. Preliminary ap- 
plications of these new atomic data and diagnostics to observations yield lower temperatures 
and higher densities than CELs, as one expects. 



1.9 Summary 

To summarize, nearly one and half centuries after William Huggins' discovery of bright 
emission lines in the spectra of gaseous nebulae, enormous progress has been achieved in 
understanding them. In particular, we believe we now understand the faint ORLs emitted by 
heavy element ions, more than half a century after Wyse's pioneer work on O II. The theory 
of photoionized gaseous nebulae seems to be solid. Yet any progress in our understanding 
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of their nature seems to be always accompanied by the emergence of new puzzles. Clearly, 
those beautiful heavenly objects, queens of the night, are more delicate than we think and 
still hold some fascinating secrets from us. With dedicated observations, more discoveries 
are sure to come. 
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